Amendments to the Specification 

Please replace paragraph [001 1] with the following amended paragraph: 

[001 1] Figs. §a-l 6a^ show a composition modulation table: Niobium 
concentration, soft local modes and mean field calculation as a result of varying 
parameters A and X, as defined hereinafter. In Figs. 6a, 6e,and 6i, modulation is by a 
sine wave with A=0.25 and A=12a. In Figs. 6b, 6f, and 6j, modulation is by a sawtooth 
wave with A=0.25 and A=12a. In Figs. 6c, and 6g , and 6k , modulation is by a sine wave 
with A=0.5 and A=12a. In Figs. 6d, and 6h , and 61 , modulation is by a sawtooth wave 
with A=0.25 and A=4a. 

Please replace paragraph [0037] with the following amended paragraph: 

[0037] Analyzing the role of the composition modulation from the 
perspective of electrostatics proves to be a valuable exercise as it leads to 
understanding the effects of different forms and parameters of the composition 
modulation. Figs. 6a-4 6aj comprise a table including: a picture of various composition 
modulation schemes, the cartesian coordinates of the local modes by plane according 
to the above given modulation, and the calculated electric field for each modulation, 
ecM(k) of Eq. 10, all indexed by plane, k. The defining parameters include the amplitude 
and the wavelength of modulation. The first modulation, as in Fig. 6a, is given by a 
sinusoidal function with an amplitude of A = 0.25, half the maximum possible fluctuation 
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amplitude. The wavelength X = 12a ensures the modulation covers twelve planes 
through one undulation. The Niobium concentration begins at 50%, becomes sparse, 
then increases and returns to 50% again, in a sinusoidal fashion. Scandium 
concentration has a similar behavior but phase shifted by 7t. The resulting local modes 
averaged over the supercell are presented below in Fig. 6e. u x , u y and u z , the average 
local modes over the supercell, are given by horizontal lines. The average local mode 
vectors are all equal, as they are in the disordered structure, and the overall electrical 
polarization then lies in the [111] direction. u x (k) and u y (k) are equal within statistical 
uncertainty, and fluctuate very little by plane. u 2 (k) sinusoidally varies around its 
average value u z as a function of k. As a result, the polarization is rotating locally with 
the fluctuation of u 2 (k). With increasing k, from 1 to 4, the polarization rotates from an 
intermediate direction between [001] and [111] towards the [1 1 1] direction. Then when 
k = 4, a value where there is no difference in the composition in the planes above, k+1 , 
and below, k-1, the polarization is that of the disordered and average structure, in the 
[111] direction. As k increases further, from 5 to 10, the polarization begins to rotate 
from the [111] direction toward the [110] direction as u 2 (k) becomes suppressed. Then 
to complete the picture, for k values 1 1 and 12, u 2 (k) begins to rise again and the 
polarization rotates through [1 1 1] towards the [001] direction. This rotation of the 
electrical polarization is quickly corroborated in the theoretical calculation of the local 
electric field, also dependent on plane index k (Fig. 6i). This is the scivi(k) given by Eq. 
8. Associated with the composition modulation is this e C M(k), an internal electric field 
capable of rotating the local polarization from that of the mean structure. More 
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precisely, a value of e C M(k) > 0 yields local mode u z (k) > u z . Likewise if e C M(k) < 0, local 
mode u 2 (k) < u z . 
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